Quantum non-demolition measurements of a qubit coupled to a harmonic oscillator 
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We theoretically describe the weak measurement of a two-level system (qubit) and quantify the 
degree to which such a qubit measurement has a quantum non-demolition (QND) character. The 
qubit is coupled to a harmonic oscillator which undergoes a projective measurement. Information 
on the qubit state is extracted from the oscillator measurement outcomes, and the QND character 
of the measurement is inferred from the result of subsequent measurements of the oscillator. We use 
the positive operator valued measure (POVM) formalism to describe the qubit measurement. Two 
mechanisms lead to deviations from a perfect QND measurement: (i) the quantum fluctuations of 
the oscillator, and (ii) quantum tunneling between the qubit states |0) and |1) during measurements. 
Our theory can be applied to QND measurements performed on superconducting qubits coupled to 
a circuit oscillator. 

PACS numbers: 03.65.Ta, 03.67.Lx, 42.50.Dv, 42.50.Pq, 85.25.-j 
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I. INTRODUCTION 

The possibility to perform repeated quantum measure- 
ments on a system with the least possible disturbance 
was first envisioned in the context of measuring gravi- 
tational wavesi. In quantum optics the optical Kerr ef- 
fect provided an early playground for studying for QND 
measurement! ] 2 ' 3 ' 4 , that were extended to the framework 
of cavity quantum electrodynamics (cavity-QED) and 
mesoscopic mechanical oscillators i 5 ' 6 ' 7 ' 8 ' 9 

The application of such a scheme to quantum infor- 
mation has stimulated great interest, in particular in 
the field of quantum computation, where fast and ef- 
ficient readout is necessary, and error correction plays 
an important role^. Schemes for qubit QND measure- 
ment have been theoretically proposed and experimen- 
tally realized with a superconducting qubit coupled to 
harmonic oscillators, either represented by external tank 
LC-circui t 11 ' 12 ! 13 : 14 ! 15 ] 16 ! 17 , or by superconducting res- 
onator that behaves as a one mode quantum harmonic os- 
cillator in circuit-QED i 18 ' 19 i 20 ' 21 i 22 i 23 i 24 A measurement 
scheme based on the Josephson bifurcation amplifier 
(JBA) 2 ^' 26 has been adopted with the aim to perform 
QND measurements of superconducting qubi t 27 ' 28 . In 
these experiments a deviation of ~ 10% from perfect 
QND behavior has been found. 

Motivated by those recent experimental achievements 
we analyze a measurement technique based on the cou- 
pling of the qubit to a driven harmonic oscillator. A 
quadrature of the harmonic oscillator is addressed via a 
projective measurement. The qubit that is coupled to 
the oscillator affects the outcomes of the measurement of 
the oscillator and information on the qubit state can be 
extracted from the results of the projective measurement 
of the oscillator. We aim to shed some light on the pos- 
sibilities to perform qubit QND measurements with such 
a setup, and to understand whether deviations from the 
expected behavior could arise from quantum tunneling 
between the qubit states. Such a tunneling process, al- 
though made small compared to the qubit energy split- 
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FIG. 1: (Color online) a) Schematics of the 4- Josephson junc- 
tion superconducting flux qubit surrounded by a SQUID, 
b) Measurement scheme: bl) two short pulses at frequency 
Ve 2 + A 2 , before and between two measurements prepare the 
qubit in a generic state. Here, e and A represent the en- 
ergy difference and the tunneling amplitude between the two 
qubit states. b2) Two pulses of amplitude / and duration 
t\ = T2 = 0.1 ns drive the harmonic oscillator to a qubit- 
dependent state, c) Perfect QND: conditional probability 
P(0j0) for A = to detect the qubit in the state "0" vs 
driving time ri and Ti, at Rabi frequency of 1 GHz. The os- 
cillator driving amplitude is chosen to be f/2n = 50 GHz and 
the damping rate k/2-k = 1 GHz. d) Conditional probability 
P(0|0) for A/e = 0.1, //2tt = 20 GHz, k/2tt = 1.5 GHz. 
A phenomenological qubit relaxation time Ti = 10 ns is as- 
sumed. 



ting, violates the QND conditions. 

One of the possible implementations of the system un- 
der consideration is the four-junction persistent current 
qubifci 4 -' 2 - 7 -' 2 ^ (flux qubit) depicted in FigQjV). It consists 
of a superconducting loop with four Josephson junctions 
and its low temperature dynamics is confined to the two 
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lowest-energy states. For an external magnetic flux close 
to a half-integer multiple of $ = h/2e, the supercon- 
ducting flux quantum, the two lowest-energy eigenstates 
are combinations of clockwise and counter clockwise cir- 
culating current states. These two states represent the 
qubit. The measurement apparatus consists of a super- 
conducting quantum interference device (SQUID), com- 
posed by two Josephson junctions, inductively coupled to 
the qubit loop. The SQUID behaves as a non-linear in- 
ductance and together with a shunt capacitance forms 
a non-linear LC-oscillator, which is externally driven. 
The two qubit states produce opposite magnetic field 
that translate into a qubit dependent effective Joseph- 
son inductance of the SQUID. The response of the driven 
SQUID is therefore qubit-dependent. 

In order to treat the problem in a fully quantum me- 
chanical way, we linearize the SQUID equation of motion, 
such that the effective coupling between the driven LC- 
oscillator and the the qubit turns out to be quadratic. 
The qubit Hamiltonian is Tis = + Acrx/2. In the 

experimental, the tunneling amplitude A between the 
two qubit current states is made small compared to the 
qubit gap E = \J t 2 + A 2 , therefore also A <C e, such 
that it can be considered as a small perturbation. The 
absence of the tunneling term would yield a perfect QND 
Hamiltonian (see below). From the experimental param- 
eters A = 5 GHz and E = 14.2 GHa2L2a, it follows that 
A/f« 0.38, yielding a reduction of the visibility in Fig.[T] 
d) on the order of 10%. 

The QND character of the qubit measurement is stud- 
ied by repeating the measurement. A perfect QND setup 
guarantees identical outcomes for the two repeated mea- 
surement with certainty. In order to fully characterize 
the properties of the measurement, we can initialize the 
qubit in the state |0), then rotate the qubit by applying 
a pulse of duration t\ before the first measurement and a 
second pulse of duration between the first and the sec- 
ond measurement. The conditional probability to detect 
the qubit in the states s and s' is expected to be indepen- 
dent of the first pulse, and to show sinusoidal oscillation 
with amplitude 1 in t%. Deviations from this expectation 
witness a deviation from a perfect QND measurement. 
The sequence of qubit pulses and oscillator driving is de- 
picted in Fig. [Ib). The conditional probability P(0|0) to 
detect the qubit in the state "0" twice in sequence is plot- 
ted versus t\ and Ti in Fig.rjJ) for A = 0, and in Fig.QJl) 
for A/e = 0.1. We anticipate here that a dependence on 
T\ is visible when the qubit undergoes a flip in the first 
rotation. Such a dependence is due to the imperfections 
of the mapping between the qubit state and the oscillator 
state, and is present also in the case A = 0. The effect of 
the non-QND term Acrx results in an overall reduction 
of P(0|0). 

In this paper we study the effect of the tunneling 
term on the quality of a QND measurement. Many at- 
tempts to understand the possible origin of the devia- 
tions from perfect QND behavior appearing in the exper- 
iments have been concerned with the interaction with the 



environmen t 22 ' 23 ^ 3 ^ 3 ^ 3 ^ 3 ^. The form of the Josephson 
non-linearity dictates the form of the coupling between 
the qubit and the oscillator, with the qubit coupled to the 
photon number operator of the driven harmonic oscilla- 
tor, aza^a, rather than to one quadrature, <Tx{a + a'), 
and the effect of the tunneling term ax present in the 
qubit Hamiltonian is considered as a small perturbation. 

The work we present is not strictly confined to the 
analysis of superconducting flux qubit measurement. 
Rather, it is applicable to a generic system of coupled 
qubit and harmonic oscillator that can find an applica- 
tion in many contexts. Moreover, the analysis we present 
is based on the general formalism of the positive opera- 
tor valued measure (POVM), that represents the most 
general tool in the study of quantum measurements. 

The paper is structured as follows: in Sec. [Til we in- 
troduce the idea of QND measurement and describe the 
conditions under which a QND measurement can be per- 
formed. In Sec. IIHI we derive the quadratic coupling be- 
tween the qubit and the oscillator and the Hamiltonian 
of the total coupled system. In Sec. IIVI we construct 
the qubit single measurement with the POVM formal- 
ism and in Sec. [V] we consider the effect of the non-QND 
term in the POVM that describes the single measure- 
ment. In Sec. I VII we construct the two- measurement 
formalism, by extending the formalism of POVM to the 
two subsequent measurement case. In Sec. IVIII we con- 
sider the single measurement in the case A = and study 
the condition for having a good QND measurement. In 
Sec. I VII J we calculate the contribution at first order and 
second order in A/e to the POVM and to the outcome 
probability for the qubit single measurement. In Sec. IIXI 
we calculate the contribution at first and second order 
in A/e to the POVM and to the outcome probability for 
the qubit two subsequent measurements. In Sec. [X] we 
study the QND character of the measurement by looking 
at the conditional probability for the outcomes of two 
subsequent measurements when we rotate the qubit be- 
fore the first measurement and between the first and the 
second measurement. 



II. QND MEASUREMENTS 

We consider a quantum system on which we want to 
measure a suitable observable A. A measurement proce- 
dure is based on coupling the system under consideration 
to a meter. The global evolution entangles the meter and 
the system, and a measurement of an observable B of the 
meter provides information on the system. In general, a 
strong projective measurement on the meter translates 
into a weak non-projective measurement on the system. 
This is because the eigenstates of the coupled system dif- 
fer in general from the product of the eigenstates of the 
measured observable on the system and those of the me- 
ter. 

Three criteria that a measurement should satisfy in 
order to be QND have been formulate d 31 ! 35 : i) correct 
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correlation between the input state and the measurement 
result; ii) the action of measuring should not alter the 
observable being measured; iii) repeated measurement 
should give the same result. These three criteria can be 
cast in a more precise way: the measured observable A 
must be an integral of motion for the coupled meter and 
system, 1 . Formally this means that the observable A that 
we want to measure must commute with the Hamiltonian 
7i, that describes the interacting system and meter, 

[H,A} = 0. (1) 

Such a requirement represents a sufficient condition in 
order that an eigenstate of the observable A, determined 
by the measurement, does not change under the global 
evolution of the coupled system and meter. As a conse- 
quence, a subsequent measurement of the same observ- 
able A provides the same outcome as the previous one 
with certainty. 

Finally, in order to obtain information on the system 
observable A by the measurement of the meter observable 
B, it is necessary that the interaction Hamiltonian does 
not commute with £?, 

[Hint, 5] ^0, (2) 

where Hint describes the interaction between the meter 
and the system, 

H = Hs + H mc tcr + Hint- (3) 

Altogether, these criteria provide an immediate way to 
determine whether a given measurement protocol can 
give rise to a QND measurement. At this level the ob- 
servables A and B and the Hamiltonian TL do not per- 
tain to any particular system. In the next section we will 
identify each term for the system we want to study. 

III. MODEL: QUADRATIC COUPLING 

As far as the application of our model to the mea- 
surement of a persistent current qubit with a SQUID is 
concerned, we provide here a derivation of the quadratic 
coupling mentioned in the introduction. 

We identify the system with a flux qubit that will be 
described by the Hamiltonian Tis- The meter is repre- 
sented by a SQUID and it is inductively coupled to a 
flux qubit via a mutual inductance, in such a way that 
the qubit affects the magnetic flux through the loop of 
the SQUID. The Hamiltonian that describes the SQUID 
and the interaction with the qubit can be written as 

H me ter + H n t = ^ - J± COS (2^$/$ ) COS (4) 

where <p — tfi — <p% is the difference of the phases of 
the two Josephson junctions ipi and i{>2 that interrupt 
the SQUID loop, Lj the Josephson inductance of the 
junctions (nominally equal), and Q is the difference of 



the charges accumulated on the capacitances C that 
shunt the junctions. Up to a constant factor, (p and Q 
are canonically conjugate variables that satisfy [<p, Q] = 
2ei. We split the external flux into a constant term 
and a qubit dependent term, such that cos(27r$/<i>o) = 
cos(27r$ oxt /$o + 2irMI q <T Z /^ ) = A + \icr z , with I q 
the current in the qubit loop and M the mutual induc- 
tance between qubit and SQUID loop. Expanding the 
potential up to second order in (/?, one obtains 

H m eter + Wi a t ~ ^ + (Ao + AlCTz) ^^-J 2L~' ^ 

with A = cos(27ri> ex t/$o) cos(27rM/ 9 /<i>o) and Ai = 
— sin(27r<I>ext/$o) sin(27rM/ 9 /<I>o). We introduce the zero 
point fluctuation amplitude a — (Lj/XqC) 1 / 4 , the bare 
harmonic oscillator frequency c^ho = i/Ao/ijC, and the 
in-phase and in-quadrature components of the field 

^<P = X = a^(a + a% (6) 

with a and harmonic oscillator annihilation and cre- 
ation operators satisfying [a, a'] = 1. Apart from a renor- 
malization of the qubit splitting, the Hamiltonian of the 
coupled qubit and linearized SQUID turns out to be 

H mo tcr + Hint = ^ho(l + ffCT Z )a ] a 

+ huj ho gcrz{a 2 + a f2 ), (8) 

with g = Ai/2A = tan(27r$ cxt /$ ) tan(27rM/ g /$ )/2. 
The frequency of the harmonic oscillator describing the 
linearized SQUID is then effectively split by the qubit. 

The Hamiltonian can now be written in the form of 
Eq. §3§ with an additional driving term (from here on we 
set h = 1), 

meter H~ TYint '^drive 

(*)■ (9) 

The qubit Hamiltonian written by means of the Pauli 
matrices <7i (we denote 2x2 matrices in qubit space with 
bold symbols) in the basis of the current states {|0), |1)} 
is 

n s = ^<Tz + j<rx, (io) 

where e = 2/ 9 (<I>ext — < &o/2) represents an energy differ- 
ence between the qubit states and A the tunneling term 
between these states. The Hamiltonian of the oscillator 
(or SQUID) is 

Hmctcr = Wharfa. (11) 

The Hamiltonian that describes the coupling between the 
qubit and the harmonic oscillator in the rotating wave 
approximation (RWA), where we neglected the terms like 
a 2 and a) 2 , is given by 

Hint = gcrza^a, (12) 
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with g = Who<7 [361 ] , and the external driving of the har- 
monic oscillator is described by 



Wdrive(t) = /(t)(a + a t ). 



(13) 



and throughout this work, we choose a harmonic driv- 
ing force f(t) = 2f cos(w^t). Neglecting the fast rotating 
terms ae _Md * and a)e luJdt , after moving in the frame ro- 
tating with frequency u>d, the Hamiltonian becomes time 
independent, 



H = Us 



•/(o + a + ), 



(14) 



with Aujz = uz-Wd, and the qubit-dependent frequency 
given by u z = w ho (l + gcr z )- 

The qubit observable that we want to measure is 
A = <tz and, due to the presence of the term Aax /2, 
it does not represent an integral of the motion for the 
qubit, [Hs,cz] ^ 0. Therefore the measurement is not 
supposed to be QND, Eq. |T]) not being satisfied. How- 
ever, for A <C e the variation in time of <Jz becomes slow 
on the time scale determined by 1/e and one expects 
small deviations from an ideal QND case. The presence 
of the non-QND term crx term in TLs inhibits an exact 
solution and a perturbative approach will be carried out 
in the small parameter A/e <C 1 . 



IV. SINGLE MEASUREMENT 

The weak measurement of the qubit is constructed as 
follows. We choose the initial density matrix (t = 0) 
of the total coupled system to be the product state 
p(0) = p ® |0)(0|, with the qubit in the unknown ini- 
tial state p and the oscillator in the vacuum state |0), 
and we let the qubit and the oscillator become entangled 
during the global time evolution. We then assume that 
at time t we perform a strong measurement of the flux 
quadrature X = <r(a + a))/^/2, by projecting the oscilla- 
tor on to the state |x)(x|. Such a state of the oscillator 
is quite unphysical, it has infinite energy and infinite in- 
determinacy of the P = (a — a^)/V2i quadrature. More 
realistically, what would happen in an experiment is that 
the oscillator is projected on to a small set of quadrature 
states centered around x. This can be described as a 
convolution of the projector |x)(x| with a distribution 
characteristic of the measurement apparatus, that can 
be included in the definition of the qubit weak measure- 
ment. However, we choose to keep the model simple and 
to work with an idealized projection. 

In the interaction picture, the projection on the state 
|x)(x| corresponds to the choice to measure the quadra- 
ture X(t) = a(ae- luJbat 



x(t) = Tr[Xp(t)} = Tr [X(t)p R (t) 
p R (t) = U R {t) p(0) 4(t). 



(15) 
(16) 



where an expression oiU R (t) and its derivation is given by 
Eq. (| A5|) in Appendix [Al The operator U R (t) describes 



the time-evolution of p in the rotating frame. The prob- 
ability to detect the outcome x can then be written as 

Prob(M) = Tr[(x\p R (t)\x)} 



= Tr 



(x\U R (t)\0) p (6\W R (t)\x) , (17) 



where the trace is over the qubit space, and {\x}} is a 
basis of eigenstates of X(t). We define the operators 



N(x,t) = (x\U R (t)\0), 
F(x,t) = N t (a;,i)N(x,t), 



(18) 
(19) 



acting on the qubit and, using the property of invariance 
of the trace under cyclic permutation, we write 



Prob(x,t) = Tr F(z,t)p(0). 



(20) 



The state of the system after the measurement is p(x, i)t 
\x)(x\, with the qubit in the state 



p(x,t) 



N(a;,t)p(0)Nt(a;,t) 
Prob(x, t) 



(21) 



The operators F(x, t) are positive, trace- and hcrmiticity- 
preserving superoperators (i.e. they map density opera- 
tors into density operators) acting on the qubit Hilbert 
space. Moreover, they satisfy the normalization condi- 
tion 



dxF(x.t) = t, 



(22) 



from which the conservation of probability follows. 
Therefore, they form a positive operator valued mea- 
sure (POVM), and we will call the operators F(x, t) a 
continuous POVM. We point out here that modeling a 
more realistic scenario, by including a convolution of the 
projector |x)(x| with a distribution characteristic of the 
measurement apparatus, corresponds to the construction 
of a more general POVM of the harmonic oscillator, that 
would not alter qualitatively the description of the qubit 
measurement in terms of POVM. 

The probability distribution Prob(ir, t) depends 
strongly on the initial qubit state p . In general 
Prob(x, t) is expected to have a two-peak shape, aris- 
ing from the two possible states of the qubit, whose rel- 
ative populations determine the relative heights of the 
two peaks, one peak corresponding to |0) and the other 
to|l). 

We now define an indirect qubit measurement that 
has two possible outcomes, corresponding to the states 
"0" and "1". As a protocol for a single-shot qubit mea- 
surement, one can measure the quadrature X and assign 
the state "0" or "1" to the qubit, according to the two 
possibilities of the outcome x to be greater or smaller 
than a certain threshold value Xth, x > Xth — * |0) or 
x < x t h — > as depicted in Fig. [21 Alternatively, we 
can infer the qubit state by repeating the procedure many 
times and constructing the statistical distribution of the 
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outcome x. We then assign the relative populations of 
the qubit states |0) and |1) by respectively integrating 
the outcome distribution in the regions 77(1 ) = [xt^oo), 
?/(-!) = (-°°>Sth)< 

We formally condensate the two procedures and define 
a two-outcome POVM, that describes the two possible 
qubit outcomes, by writing 

F(s,i) = / dxF(x,t), (23) 

J V (s) 

Prob(s,i) = Tr[F(s,i)p(0)], (24) 

with s = ±1. We will call F(s,t) a discrete POVM, in 
contrast to the continuous POVM F(a;, t) defined above. 
Here, we introduce a convention that assigns s = +1 to 
the "0" qubit state and s = — 1 to the "1" qubit state. 
The probabilities Prob(s, t) are therefore obtained by in- 
tegration of Prob(x,t) on the subsets rj(s), Prob(s,i) = 
J^/g) dxProb(x, t). On the other hand, the probability 
distribution Prob(ir, t) is normalized on the whole space 
of outcomes which leads to P(0,t) + P(l,t) — 1 at all 
times. Typically, it is not possible to have a perfect map- 
ping of the qubit state. 

V. EFFECTS OF THE TUNNELING a x TERM 

Deviations from an ideal QND measurement can arise 
due to the presence of a non-zero ax term in the qubit 
Hamiltonian. In SC flux qubits, such a term is usually 
present; it represents the amplitude for tunneling through 
the barrier that separates the two wells of minimum po- 
tential, where the lowest energy qubit current states are 
located. This term cannot be switched off easily. 

We can expand the full evolution operator Un(t) in 
powers of At, as in Eq. (|A11[) , and obtain a formally 
exact expansion of F(x,t), 

00 

F(i,i) = ^FW(M). (25) 

n=0 

Due to the transverse (X _L Z) character of the pertur- 
bation it follows that the even terms in this series (cor- 
responding to even powers of At) have zero off-diagonal 
entries, whereas the odd terms have zero diagonal entries. 
Due to the normalization condition Eq. (f2"2")l , valid at all 
orders in At, it can be shown that 

J dxF^(x, t) = S n , l, (26) 
and consequently, 

F^(s,t) = 6 n , l. (27) 

s=±l 

As a result, the probability Prob(s,i) is given as a power 
expansion in the perturbation 

00 

Prob(s, t) = ^ Prob (n) (s, t), (28) 

n=0 



where £ s =±i Prob (n) (s, t) = <$„,„. 

The expansion of the evolution operator and conse- 
quently of the continuous and discrete POVMs is in the 
parameter At. The requirement that the deviations in- 
troduced by the tunneling <r x term in the time-evolution 
behave as perturbative corrections sets a time scale for 
the validity of the approximation, namely t <C 1/A, for 
which we will truncate the expansion up to second or- 
der. The tunneling non-QND term is considered as a 
perturbation in that experimentally one has A/e -C 1. 
It turns out to be convenient to choose as a time scale 
for the qubit measurement t ~ 1/e, for which follows 
At ~ A/e < 1. 

VI. TWO SUBSEQUENT MEASUREMENTS 

A QND measurement implies that repeated measure- 
ments give the same result with certainty. In order to 
verify such a property of the measurement, we construct 
here the formalism that will allow us to study the corre- 
lations between subsequent measurements. 

After the oscillator quadrature is measured in the first 
step at time t and the quadrature value x is detected, the 
total system composed of the qubit and the oscillator is 
left in the state p(x,t) ® The fact that we can 

split the total state after the measurement into a prod- 
uct state is a consequence of the assumption that the 
measurement of the harmonic oscillator is a projection. 
Had a more general POVM of the harmonic oscillator be 
involved, then such a conclusion would not hold. After 
the first measurement is performed, the total system is 
left alone under the effects of dissipation affecting the os- 
cillator. A harmonic oscillator which is initially prepared 
in a coherent state evolves, under weak coupling to a 
bath of harmonic oscillators in thermal equilibrium, to a 
mixture of coherent states with a Gaussian distribution 
centered around the vacuum state (zero amplitude co- 
herent state) with variance nth = (exp(fiw/fcsT) — l) -1 , 
uj being the frequency of the harmonic oscillator, T the 
temperature, and ks the Boltzmann constant, whereas 
in the case T = it evolves coherently to the vacuum |0) 

We now assume that the state of the total system 
(qubit and oscillator) before the second measurement is 

p(x,*)(g)|6)<0|. (29) 

Following the previously described procedure for the 
qubit single- measurement, a second measurement of the 
quadrature X with outcome y performed at time t', hav- 
ing detected x at time t, would yield the conditional prob- 
ability distribution 

Prob(y, t'\x, t) = Tr [F(y, t')p(x, t)} . (30) 

Defining the continuous POVM qubit operators for two 
measurements as 

F(y,t';x,i) = -N\x,t)F(y,t' -t)N(x,t), (31) 
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FIG. 2: (Color online) Schematic description of the sin- 
gle measurement procedure. In the bottom panel the 
coherent states |ao) and \oci), associated with the qubit 
states |0) and |1), are represented for illustrative pur- 
poses by a contour line in the phase space at HWHM 
of their Wigner distributions, defined^ as W(a,a*) — 
(2/7r 2 )exp(2|a| 2 ) J d/3{— j3\p\/3) exp(/3a* — f3*a). The corre- 
sponding Gaussian probability distributions of width a cen- 
tered about the qubit-dependent "position" x s are shown in 
the top panel. 



the joint probability distributions for two subsequent 
measurements is 

Prob(y, if; x, t) = Prob(y, t'\x, t)Prob(a;, t) (32) 
= Tr[F(y,t';x,t)p }. (33) 

The operators F(y, t';x, i) satisfy the normaliza- 
tion condition J dx J dyF(y,t';x,t) = 1, ensur- 
ing the normalization of the probability distribution 
J dx J dyProb(y,t';x,t) — 1. By inspection of Eqs. 
and A3 ID. it follows that 



J dyF(y,t';x,t)=F(x,t), 



(34) 



and the marginal distribution for the first measurement 
is 

Prob M (M) = / dy Pioh(y,t';x,t) = Tr[F{x, t)p ], 

(35) 

stating that the probability to detect x in the first mea- 
surement is independent on whatever could be detected 
in the second measurement. On the other hand, the 
marginal probability distribution for the second measure- 
ment turns out to be 

Prob M (y,t') = J dxProb(y,t';x,t)=Tr[F{y,t'-t)p(t)}, 

(36) 



where p(t) = Tr s [U R (t)p ® \6)(6\U R (t)} is the qubit re- 
duced density matrix at time t. We define the discrete 
POVM for the correlated outcome measurements as 

F(s',t';s,t)= [ dx f dyF(y,t';x,t). (37) 

Analogously to Eq. ([34]) it follows that F(s,t) = 
Y^ s i F(s', t'; s, t), and the probability distribution for the 
outcomes of the two subsequent measurement is simply 
given by 



Prob(s', f'; s, t) = Tr[F(s', t'\ s, t)p ] 



(38) 



and it follows that Prob(s', t'\ s, t) — Prob(s, t) — 
Tr[F(s, t)p ). The conditional probability to obtain a 
certain outcome s' at time t', having obtained s at time 
t, is given by 

, , , „ , Tr\F(s't';s,t)p () \ , s 

The discrete POVM for the double measurement can 
be in general written as 

F{s',t'; s, t)=- [F(s', t')F(s, t) + h.c] + C(s', t'- s, i), 

where we have symmetrized the product of the two 
single-measurement discrete POVM operators F(s',t') 
and F(s, t) in order to preserve the hermiticity of each of 
the two terms of Eq. (j40|) . 

Proceeding as for the case of a single qubit measure- 
ment, we expand F(y, t'; x, t) in powers of A/e. Equating 
all the equal powers of A/e in the expansion it follows 
that 

F(")(M) = E F( "V,i';M), (41) 

s' 

with E s yF (n) (s'i';M)=<W)l. 

VII. IDEAL SINGLE MEASUREMENT 

The dynamics governed by (t) produces a coher- 
ent state of the oscillator, whose amplitude depends on 
the qubit state, see Fig. O In this case the continu- 
ous POVM operators have the simple form F^(x,t) = 
(atz(t)\x}(x\az(t)}, defined through Eq. (|A9[) in the Ap- 
pendix lAl In the er^-diagonal basis {\i}}, with i = 0, 1, 
it is given by 



F^ix^i^SijGix-Xiit)), 



(42) 



where Xi{t) = y/2aRje[ai(t)] and Q{x) = 
cxp(— x 2 /a 2 )/<Ty/Tr is a Gaussian of width a schemat- 
ically depicted in Fig. O Introducing a rate k that 
describes the Markovian damping of the harmonic oscil- 
lator by a zero-temperature bath of harmonic oscillators, 
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the coherent state qubit-dependent amplitude aii(t) is 
found to be2£ 



Oi(t) = Ate*** 



1 



-iAuit-Kt/2 



(43) 



with Auji — LL>i — LUd and the qubit-dependent amplitudes 
and phases given by 



Ai = 



f 



« 2 /4 



arctan 



k/2 



7T 

2' 



(44) 
(45) 



The probability distribution for the X quadrature 
outcomes is then given by the sum of the two qubit- 
dependent Gaussians, weighted by the initial state state 
occupancy, and the discrete POVM for the qubit mea- 
surement as given by Eq. (|23p becomes 



F<°)(M) 



1 + s erf 



^ 6x(t) 



crz 



(46) 



where s = ±1 labels the two possible measurement 
outcomes, and Sx(t) = cRe 6a(t) /y/2, where Sa(t) = 
ao(t) — a\{t). The indirect qubit measurement gives the 
outcome probability 



Prob( S ,i) = - 



erf 



8x(t) 



{crzja 



(47) 



with ((Tz)o = Tt[ctzPq]- Supposing that the qubit is 
prepared in the |0) state, one expects to find Prob(O) = 1 
and Prob(l) = 0. From Eq. (|47[) . we see that even for 
A = this is not always the case. 



A. Short time 

We choose a time t « 1/e and a driving frequency close 
to the bare harmonic oscillator frequency. We can then 
expand the qubit dependent signal and obtain the short 
time behavior of the signal difference, 



Sa(t) ^ V2tA, 



(48) 



with A = f(e 2i <t>° - e 2 *)/\/2. The first non-zero con- 
tribution is linear in t, because the signal is due to the 
time-dependent driving™ We measure a rotated quadra- 
ture X v = a(ae^ tv + a^e lip )/V2, and choose the phase 
of the local oscillator such that <p = arg A. With this 
choice we have 8x(t) — o\A\ t, and the probabilities for 
the two measurement outcomes 



Prob(M) = i[l 



8(tr z ) Q ed(\A\ t)] 



(49) 



In Fig. [3] we plot the probability of measuring the "0" 
state Prob(0,i = 0.1 ns) as a function of the detuning 
Auj = Who — LOd and the driving amplitude /, given that 



X 

K 

& 




Aco/2jr (GHz) 

FIG. 3: (Color online) Prob(0, t — 0.1 ns) for the initial state 
|0}(0|, as given by Eq. (|49|l . plotted as a function of the de- 
tuning Auj/2n. The values of the parameters used are listed 
in Tab. U 



the initial state is "0", p = |0)(0|. It is possible to 
identify a region of values of / and Aw where erf(|A|t) ~ 
1 HH. It then follows that 



Prob(s) 



1 



s{(tz)o\ 



(50) 



This case corresponds to a strong projective measure- 
ment, for which the outcome probabilities are either 
ore 1, thus realizing a good qubit single measurement. 

For driving at resonance with the bare harmonic os- 
cillator frequency cuho, the state of the qubit is encoded 
in the phase of the signal, with <j>\ = —cj)Q, and the am- 
plitude of the signal is actually reduced, as also shown 
in Fig. [3] for Aw = 0. When matching one of the two 
frequencies u>i the qubit state is encoded in the ampli- 
tude of the signal, as also clearly shown in Fig. [3] for 
Auj = ±g. Driving away from resonance can give rise to 
significant deviation from and 1 to the outcome prob- 
ability, therefore resulting in an imprecise mapping be- 
tween qubit state and measurement outcomes and a weak 
qubit measurement. 



VIII. CORRECTIONS DUE TO TUNNELING 

In order to compute the correction at first order in the 
tunneling term proportional to A we expand the evolu- 
tion operator Un(t) up to first order in At. By making 
use of the expression Eq. (|A12|) for the perturbation in 
the interaction picture, the off-diagonal element of the 
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FIG. 5: (Color online) Plot of the second order correction 
Prob (2) (s = 1) to detect "1" for the initial state |0> <0| , for 
At = A/e , as a function of the detuning Auj/2n, for several 
values of the driving amplitude /. The values of the parame- 
ters for the evaluation used are listed in Tab. [I] 



FIG. 4: (Color online) Probability to detect the outcome s = 
1, corrected by a) the real part of for the initial state 

|+)x(+|i and b) the imaginary part of F^ 1 ' , for the initial 
state \+)y{+\, plotted versus the detuning Alu/2tt for several 
values of the amplitude /. The values of the parameters used 
arc listed in Tab. HI 



first order correction to F(x,t) is given by 

F^(x,t) 01 = f dt' [Q (x - x (t) + 6z(t')) 

* Jo 

- g(x- Xl (t)-6z(t')*)]e iet 'r(t'), (51) 

with the complex displacement Sz(t) — aSa(t)/^/2 and 
the overlap T(t) = (a (t)\ai(t)} , where 

r(i) = exp ^\Sa(t)\ 2 - iIm[aS(i)oi(*)]) ■ (52) 

Here the state "0" is labeled by its crz- eigenvalue s = 1, 
whereas the state "1" by its az~ eigenvalue s = — 1. Anal- 
ogously to the unperturbed case, the first order contribu- 
tion to the discrete POVM is obtained by integrating the 
continuous POVM in x over the subsets i](s). Defining 
the function 

^)( f ) = 4jrV e ^r(Oerf( fa ^; fe ^), (53) 

we can write the first order contribution to the discrete 
POVM as 

F (1) (s,i) = s (Re F^it) (T X - Im F (1 \t) cry) , (54) 

and the resulting first order correction to the probability 
follows directly from Eq. ([217)1 . This correction is valid 



only for short time, t <^ 1/A. For times comparable 
with 1/A a perturbative expansion of the time evolu- 
tion operator is not valid. Choosing t w 1/e, we can 
effectively approximate the phase associated with two 
different coherent states as lm[ao(t)ai(t)*} ss ipt 2 , with 
ip = / 2 sin(20 o -200, the expression for F^(t) further 
simplifies, 

FW(t) = ^^ f dt'e^'-i^-^exi QA\(t - t')) . 

(55) 

We study the behavior of F^\t) as a correction 
to a qubit projective measurement, that is in the 
range of driving amplitudes and frequencies that ensure 
erf(|A| t) « 1. 

The real and imaginary part of F^'^t) represent the 
first order correction to the outcome probability of the 
measurement for two particular initial states, respec- 
tively and |+)y(+|, with \±) x = (|0) ± \l))/y/2 

and |±)y = (|0) ± i|l))/V2. In the first case we have 

Prob(s, t) = -+s ReF {1) (i), (56) 

and analogously for the second case, with the imaginary 
part instead of the real one. We see that the probability 
to obtain "0" is increased by ReFW(t) and the probabil- 
ity to obtain " 1" is decreased by the same amount. Since 
the contribution to first order in At only affects the off- 
diagonal elements of p , there is no effect, at first order 
for the qubit basis states |0) and |1). 

In Fig. |4] a) we plot the probability to detect the out- 
come state "0", corresponding to the outcome s = 1, 
corrected up to first order in the perturbation for At — 



9 



A/e = 0.1, for the initial state p = |+)x(+|, that in- 
volves Re FW(t). We sec that the effect of the tunnel- 
ing is largest when driving at resonance with the two 
qubit-shifted frequencies, Alu w ±g. For weak driving 
amplitude /, the phase ip in Eq. (|5 5 j) is small and the 
response is of order of ~ 1%, close to the qubit-split fre- 
quency. By increasing the strength of the driving we see 
that the structure acquires two local minima in proxim- 
ity of the resonance Aw w ±g and a maximum exactly 
at resonance Aw = ±g. The strong oscillatory behav- 
ior of the probability is due to a rapid change of sign of 
the phase ip in proximity of the qubit-split frequencies, 
that is enhanced when the driving strength / increases. 
In Fig. 2] b) we plot the probability to detect the out- 
come state " 0" , corresponding to the outcome s = 1 for 
the initial state p Q = \+)y(+\, that involves Im F^'^t). 
In comparison to Fig. [4] a) we find twice as many oscil- 
lations in the structure, typical for the imaginary part 
of a response function, when compared to the real part, 
and an overall scale factor of order 0.1. Besides, the sign 
of the response is not unique. The scale factor and the 
sign are understood from the Hamiltonian of the qubit, 
H = (ecrz + Acrx)/2, with the condition A < t. Un- 
der free evolution for a time t ~ 1/e, the initial state 
|+)x{+| acquires a larger component in the Z-direction 
than the initial state |+)y(+|. Deviation from this naive 
picture due to the coupling with the measurement ap- 
paratus translates in fluctuations that may determine a 
change of sign in the response for the case p Q = \+)y(+\- 
Away from the resonances we see no significant contribu- 
tion to the outcome probability. 

First order effects in the tunneling cannot be respon- 
sible for qubit flip during the measurement. In order to 
estimate the deviation from a perfect QND measurement 
for the eigenstates of crz, we have to consider the effect 
of the perturbation at second order. We define F^ (t) in 
Eq. (|B5[) and the contribution at second order in At to 
the discrete POVM is then 



F^(s,t) =-sF^(t) 



crz- 



(57) 



The dependence on s factorizes, as expected from the 
symmetry between the states |0) and |1), in the picture 
we consider with no relaxation mechanism. The correc- 
tion at second order in A/e to the outcomes probability 
is given by 

Prob< 2 > ( S ,t) = -sFW(t){cr z ) a . (58) 
In Fig. [5] we plot the second order correction to the prob- 



Quantity 


Symbol 


Value for plots 


Qubit detuning 


e 


2tt x 10 GHz 


Damping rate 


K 


2tt x 0.1 GHz 


Coupling strength 


9 


2tt x 0.3 GHz 


Qubit tunneling 


A/e 


0.1 



ability to obtain "1" having prepared the qubit in the 
initial state p = |0)(0|, corresponding to F^(t), for 
At = A/e = 0.1. We choose to plot only the devia- 
tion from the unperturbed probability because we want 
to highlight the contribution to spin-flip purely due to 
tunneling in the qubit Hamiltonian. In fact most of 
the contribution to spin-flip arises from the unperturbed 
probability, as it is clear from Fig. [3) Around the two 
qubit-shifted frequencies, the probability has a two-peak 
structure whose characteristics come entirely from the 
behavior of the phase ip around the resonances Alu rs ±g. 
We note that the tunneling term can be responsible for 
a probability correction up to ~ 4% around the qubit- 
shifted frequency. 

From the analysis of the qubit single measurement in 
Fig. [3] we conclude that a weak POVM qubit measure- 
ment, that yields a large error in the determination of the 
qubit state, can arise when weakly driving the harmonic 
oscillator. Therefore, only a strong qubit projective mea- 
surement, obtained for strong driving of the oscillator, 
can produce a confident mapping of the qubit state at 
the level of a single measurement. In this case, a devia- 
tion on the order of a few percent in the state assignment 
can be ascribable to the tunneling term. 



IX. QND CHARACTER OF THE QUBIT 
MEASUREMENT 



As explained in Sec. [TT1 repeated measurements should 
give the same result if the measurement is QND. Such 
a requirement means that if a measurement projects the 
system onto an eigenstate of the measured observable, 
then a subsequent measurement should give the same re- 
sult with certainty. The presence of a term that does not 
satisfy the QND condition may affect the character of 
the measurement essentially in two ways: i) by introduc- 
ing deviations from the projection character of the single 
measurement, and ii) by generating non-zero commuta- 
tors in the two-measurement POVM. These may affect 
the two-outcome probabilities. 



A. A = case 

The case A = satisfies the requirement for a QND 
measurement of the qubit observable crz- The discrete 
POVM factorizes in this particular case, by virtue of the 
fact that [N(°)(j/,t' -t),N^(x,t)} =0, 



F<°> (s', t'; s, t) = F^ (a', t' - i)F<°> (s, t). 



(59) 



Choosing, e.g., t' = It and using Eq. ([38]). the joint prob- 
ability for the two measurements reads 



Prob(s';s) 



1 + s's erf 



(5x{t) 



TABLE I: Values of the parameters used in the plots. 



V a 

+ ( S ' + S )erf(^)( 



c z)o 



(60) 
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In the region of driving frequency and amplitude that 
ensure erf (Sx/a) ~ 1, we find 

r> wii l + s's + (s + s')(a z )o ra1 s 

Prob (a s) = ~ . r . , 61 

2(1 + s{er z )o) 

with Prob(s; s) = Prob(s), and Prob(— s; s) = 0, 
and the conditional probability is Prob(s|s) = 1, and 
Prob(— s\s) — 0, regardless of (crz)o- However, it has to 
be noticed that in the case the condition eri(5x/a) w 1 
does not perfectly hold, the conditional probability for 
the two measurement to give the same outcome becomes 

2(1 + s ert{bx/a){(Tz)o) 
and this does depend on the initial state {cr z)a- 



B. First order contribution 

We now apply the perturbative approach in At to esti- 
mate the effect of the non-QND term for the joint and the 
conditional probabilities. Due to the transverse nature of 
the perturbation, it is possible to show that all the odd 
terms have off-diagonal entries, whereas even ones are di- 
agonal. At first order in At the off-diagonal term of the 
discrete POVM is given by 

F« (s',t';s,t) = |f« (t) + S -F^ (f - t ) + S 'C« (f ; t) , 

(63) 

with the quantity C^(t';t) given by Eq. (|P36|) in Ap- 
pendix [B] For the particular choice t' = 2t, for which 
the two measurement procedures are exactly the same, 
the joint probability for the initial state p Q = \+)x(+\ is 
given at first order in At by 

Prob(s',s) = i^l + ^serf^^H 

+ l(s + s')Re F« (t) + s' Re (2t; t) 

(64) 

We immediately observe that the probability is not sym- 
metric with respect to s and s' . Although the driv- 
ing times are the same, something is different between 
the first and the second measurement, and the proba- 
bility to obtain different outcomes s' = — s is different 
from zero. An analogous result holds for the initial state 
p = |+)y(+|i with the imaginary part instead of the 
real one. Now, no matter the sign of C' 1 ', the prod- 
uct — s CW is negative in one case (s = ±1). In order 
to ensure that probabilities are non-negative one has to 
choose At small enough such that the first order negative 
corrections due to remains smaller than the unper- 
turbed probability. If At is too large, one needs to take 
higher orders into account which should then ensure an 
overall non- negative probability. The behavior of 



as a function of the detuning Acu and the driving am- 
plitude / is very similar to that of and we choose 
not to display it. The only main difference arises in the 
magnitude, for which we have |CW| < \F^\. It is the 
clear that the main deviations in the two measurement 
probabilities are mainly due to the errors in the first or 
second measurement. 



C. Second order contribution 

The contribution to the discrete POVM at second or- 
der in At can be divided into a term that factorizes the 
contributions of the first and the second measurements, 
as well as a term that contains all the non-zero commu- 
tators produced in the rearrangement, 

F^( S ',t';s,t) = F(°)( S ,t)F( 2 V,<'-t) 
+ F^(M)F (0) (s',<'-i) 
1 



F( 1 )(s,i)F( 1 )(s , ,t' -t) + h.c. 



+ C( 2 >( S ',i'; S ,i). 



(65) 



The full expression of the at second order is rather 
involved. Choosing t' = 2t we then obtain 



C^(p',2t;p,t) s 



p ps 



C^(t) 



pp 



F (1) (t) , (66) 



with C (2) (£) given by Eq. jB7| in AppendixEJ The prob- 
ability to obtain identical outcomes does depend on the 
outcome s itself, and this reflects the fact that the joint 
probability still depends on the initial states of the qubit. 
At the same time, the probability for obtaining different 
outcomes does not depend on s, as expected. However, 
direct evaluation of the function C^(p', 2t;p, t) shows 
that its contribution to the probability is of order 0.1% 
and can be neglected. 



X. RABI OSCILLATIONS BETWEEN 
MEASUREMENTS 

In order to gain a full insight in the QND character of 
the measurement, we analyze the behavior of the condi- 
tional probability to detect the outcomes s and s' in two 
subsequent measurements when we perform a rotation of 
the qubit between the two measurements. Such a pro- 
cedure has been experimentally adopted in the work of 
Lupa§cu et al. [27]. When changing the qubit state be- 
tween the two measurements, only partial QND behavior 
is expected. In addition to this, we apply an initial ro- 
tation to the qubit, such that a wide spectrum of initial 
states is tested. Ideally, the complete response of this 
procedure is supposed to be independent on the time T\ , 
during which we rotate the qubit before the first measure- 
ment, and to depend only on the time T2, during which wc 
rotate the qubit between the first and the second mea- 
surements, with probabilities ranging from zero to one 
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FIG. 6: (Color online) Conditional probability to obtain a) 
s' = s = 1, b) s' = —s — 1, c) s' = — s = — 1, and d) 
s' = s = -1 for the case At = A/e = 0.1 and T\ = 10 ns, 
when rotating the qubit around the y axis before the first 
measurement for a time n and between the first and the sec- 
ond measurement for a time ti , starting with the qubit in the 
state |0}(0|. Correction in At are up to second order. The 
harmonic oscillator is driven at resonance with the bare har- 
monic frequency and a strong driving together with a strong 
damping of the oscillator are assumed, with f/2n = 20 GHz 
and k/2tx = 1.5 GHz. 



as a function of t^. Such a prediction, once confirmed, 
would guarantee a full QND character of the measure- 
ment. 

In Fig. [T]c) we plot the conditional probability -P(0|0) 
for the case A = 0, when strongly driving the harmonic 
oscillator at resonance with the bare harmonic frequency, 
Au = 0. The initial qubit state is chosen to be |0)(0|. 
No dependence on t± appears and the outcomes s and 
s' play a symmetric role. This is indeed what we expect 
from a perfect QND measurement. In Fig. [5] we plot the 
four combinations of conditional probability P(s'\s) up to 
second order corrections in At = A/e = 0.1 and with a 
phenomenological qubit relaxation time T\ — 10 ns. We 
choose Auj = 0, that is at resonance with the bare har- 
monic frequency. The initial qubit state is |0)(0|. Three 
features appear: i) a global reduction of the visibility of 
the oscillations, ii) a strong dependence on t\ when the 
qubit is completely flipped in the first rotation and iii) 
an asymmetry under change of the outcome of the first 
measurement, with an enhanced reduction of the visibil- 
ity when the first measurement produces a result that 
is opposite with respect to the initial qubit preparation 
|0)(0|. Furthermore, we find a weak dependence of the 
visibility on t\. 



FIG. 7: (Color online) Comparison of the deviations from 
QND behavior originating from different mechanisms. Con- 
ditional probability P(0|0) versus qubit driving time n and 
T2 starting with the qubit in the state |0)(0|, for a) A = 
and Ti = oo, b) A = and T\ = 2 ns, and c) A = 0.1 e and 
Ti = co. The oscillator driving amplitude is //2-7T = 20 GHz 
and a damping rate k/2tv = 1.5 GHz is assumed. 



We now investigate whether it is possible to identify 
the contribution of different mechanisms that generate 
deviations from a perfect QND measurement. In Fig. 
we study separately the effect of qubit relaxation and 
qubit tunneling on the conditional probability P(0|0). In 
Fig- 0a) we set A = and T\ — oo. The main feature ap- 
pearing is a sudden change of the conditional probability 
P — > 1 — P when the qubit is flipped in the first rotation. 
This is due to imperfection in the mapping between the 
qubit state and the state of the harmonic oscillator, al- 
ready at the level of a single measurement. The inclusion 
of a phenomenological qubit relaxation time T\ = 2 ns, 
intentionally chosen very short, yields a strong damping 
of the oscillation along t-i and washes out the response 
change when the qubit is flipped during the first rotation. 
This is shown in Fig. b) . The manifestation of the non- 
QND term comes as a global reduction of the visibility 
of the oscillations, as clearly shown in Fig. [7]c). 

At this level it is clearly possible to associate the ob- 
served features to different originating mechanisms: i) 
qubit tunneling yields an overall reduction of the visibil- 
ity of the oscillations and an asymmetry under change 
of the outcome of the first measurement qubit relaxation 
and deviation from projective measurements, ii) qubit 
relaxation results in damping along T2 and weak depen- 
dence of the oscillations on n, and iii) deviations from 
projective measurements show up mostly when the qubit 
is flipped during the first rotation. 
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The combined effect of the quantum fluctuations of 
the oscillator together with the tunneling between the 
qubit states is therefore responsible for deviation from a 
perfect QND behavior, although a major role is played, 
as expected, by the non-QND tunneling term. Such a 
conclusion pertains to a model in which the qubit QND 
measurement is studied in the regime of strong projec- 
tive qubit measurement and qubit relaxation is taken into 
account only phenomenologically. We compared the con- 
ditional probabilitiesplotted in Fig. [S] and Fig. [7] directly 
to Fig. 4 in Ref. [27l| . where the data are corrected by 
taking into account qubit relaxation, and find good qual- 
itative agreement. 

Our findings can also be compared to the experiment 
(28|, in which the QND character of the measurement is 
addressed by studying a series of two subsequent mea- 
surements, but no qubit rotation is performed between 
the two measurements. The data in Ref. [28| are affected 
by strong qubit relaxation. However, from the analy- 
sis of the joint probabilities of the outcomes of the two 
measurements provided in Ref. (28l |. one can extract the 
conditional probabilities P(0|0) ~ 83% (when starting 
with the qubit initially in the ground state and compara- 
ble to Fig. E] a) at n = r 2 = 0), and P(0|0) - 77% (after 
a 7r-pulse is applied to the qubit initially in the ground 
state, that is comparable to Fig. [6] a) at t\ — 0.5 ns and 
t% = 0). In these cases one would expect a conditional 
probability of order 1 and a weak dependence on qubit 
relaxation. A deviation of order ~ 20% can be under- 
stood within the framework of our model as arising from 
the non-QND term and from a weak qubit measurement. 
Besides, from the data provided in Ref. [28|, one can ex- 
tract a probability of ~ 17% to obtain the excited state, 
when starting with the qubit in the ground state, already 
at the level of the single measurement. Such a behavior 
cannot be understood as a result of qubit relaxation and 
it can be ascribed to deviations from a projective qubit 
measurement. 



XI. CONCLUSION 

In this paper we have analyzed the QND character of a 
qubit measurement based on coupling to a harmonic os- 
cillator that works as a pointer to the qubit states. The 
Hamiltonian that describes the interaction between the 
qubit and the oscillator does not commute with the qubit 
Hamiltonian. This would in principle inhibit a QND mea- 
surement of the qubit. The term in the qubit Hamilto- 
nian that gives rise to the non-zero commutator is small 
compared with the qubit energy gap, and in the short 
time qubit dynamics it can be viewed as a small per- 
turbation. The perturbative analysis carried out for fast 
measurements leads us to the conclusion that the effect 
of the non-QND term can manifest itself as a non neg- 
ligible correction. A perfect QND measurement guar- 
antees perfect correlations in the outcomes of two sub- 
sequent measurements, therefore QND character of the 



measurement is understood in terms of deviations from 
the expected behavior. Corrections to the outcome prob- 
abilities have been calculated up to second order in the 
perturbing term. 

The ground and excited states of the qubit are affected 
only at second order by the perturbation, but a general 
measurement protocol should prescind from the state be- 
ing measured. Therefore, in the spirit of the experiment 
of Lupascu et al. [27} , we have studied the conditional 
probability for the outcomes of two subsequent measure- 
ments when rotating the qubit before the first measure- 
ment and between the first and the second measurement. 
In the case where the QND condition is perfectly satis- 
fied, that is when the perturbation is switched off, we 
see no dependence of the conditional probability on the 
duration of the first rotation appears and the Rabi os- 
cillations between the two measurement range from zero 
to one. This behavior shows perfect QND character of 
the qubit measurement. On the other hand, the main ef- 
fects of the non-QND term manifests itself as an overall 
reduction of the visibility of the oscillations and a asym- 
metry between the outcomes of the measurements. An 
additional dependence on the duration of the first qubit 
rotation may appear if a projective measurement of the 
qubit is not achieved already in absence of the perturbing 
non-QND term. Experimentally the measurement is not 
projective and relaxation processes inhibit a perfect flip 
of the qubit before the first measurement. 

We point out that our analysis is valid only when the 
non-QND term Aerx can be viewed as a perturbation, 
that is for short time At < 1 and when the qubit dy- 
namics is dominated by the term ecrz, for A/e <C 1. Our 
analysis is not valid for the case e = 0. In the present 
study we have neglected the non-linear character of the 
SQUID, which is not relevant to the fundamental issue 
described here, but plays an important role in some mea- 
surement procedure o 25 i 26 ' 27 i 28 . 

A way to improve the QND efficiency would be simply 
to switch the tunneling off. In the case of superconduct- 
ing flux qubit, a possibility toward smaller A could be to 
gate the superconducting islands between the junctions 
of the qubit loop [40| |. As an operational scheme one 
could think of working at finite A for logical operations 
and then at A = for the measurement. 

We acknowledge funding from the DFG within SPP 
1285 " Spintronics" and from the Swiss SNF via grant 
no. PP02-106310. 



APPENDIX A: EXACTLY SOLVABLE CASE: 

A = 

In order to determine the evolution governed by the 
Hamiltonian Eq. (fl"4|) we single out the term Hq diagonal 
in the {|s, n)} basis, with \s) the eigenstates of crz and 
\n) the oscillator Fock states, 

H = H + f{a + a^) + ^<T X , (Al) 
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with Ho = ecz/2 + Auza^a. We then work in the in- 
teraction picture with respect to Hq. The Heisenberg 
equation for the density operator reads pj = —i [Hi, pj], 
with 

Hi = Hf ] +Vi, (A2) 
Hf ] = f(ae~ lA " zt +a)e lA " zt ), (A3) 

Vi = y (e iA "V+ + e- ifi »V_) , (A4) 

where we define fi n = e+2<7a T a, and <t± = {<Jx ±io*y)/2. 
We will call the evolution operator generated by 

Wj. 

The evolution operator is given by hi (t) = 
exp(— iujdta^a — iHot)hii(t). For the measurement pro- 
cedure so far defined we are interested in the evolution 
operator in the frame rotation at the bare harmonic os- 
cillator frequency. Therefore 

U R {t) = exp(-ieta z /2 - M in tt)Ui(t). (A5) 

For the case A = the model is exactly solvable 
and uf\t) can be computed via a generalization of the 
Baker-Hausdorff formula^, 

uf\t)=D{ lz {t)), (A6) 

with the qubit-dependet amplitude r Y z (t) = 
—if f Q dse lAu>zS . The operator D(a) = exp(a t a - aa*) 
is a displacement operator—, and it generates 
a coherent state when applied to the vacuum 
| a) = D(a)\0) = e-l a l 2 / 2 X) n (Q! n /Vn!)|n). In the 
frame rotating at the bare harmonic oscillator frequency, 
the state of the oscillator is a coherent state whose 
amplitude depends on the qubit state. A general initial 
state 

Ptot(0)= E Ptfl*>0'l®|6)(6|, (A7) 

y=0,l 

where |0) is the harmonic oscillator vacuum state, evolves 
to 

P R (t)= E Pij\i){j\®\<t)){<Xj(t)l (A8) 

ii=0,l 

where we define the qubit operators ot z (t) = 
j z (t)e~ igtrTz , and the object 

\a z (t)) = D(a z )\b), (A9) 

that gives a qubit-dependent coherent state of the har- 
monic oscillator, once the expectation value on a qubit 
state is taken, = (i\a z (t)\i), for % = 0, 1. 

1. Perturbation theory in A 

For non-zero A, a formally exact solution can be writ- 
ten as 

Ui(t)=uP(t)Texp (-i A dt'Vi(t') \ , (A10) 



with Vj(t) = Uf^{t)Vi{t)Uf'{t) and T the time order 
operator. For a time scale t <C 1/A we expand the evo- 
lution operator in powers of At <C 1 , 

Ui{t) w hlf\t)(t-iAt J dsVi{st) 

- {At) 2 J ds J ds'Vi(s t)Vi(s' t)\ .(All) 
The interaction picture potential can be written as 

Vi{t)= 1 -[V{t)a + + V\t)<r-], (A12) 
with the oscillator operators T>(t) defined as 

v{t) = z>t (7o(t))e «„ tjD(7i(i)) (A13) 

= exp (iet — Hm[ao(t)ai(t)*]) 

x £>(-<5a(i)e 49 *)e 249 * Q+a . (A14) 

Here 5a(t) — a (t) — a\{t) is the difference between the 
amplitudes of the coherent states associated with the two 
possible qubit states. 



APPENDIX B: FIRST AND SECOND ORDER 
QUANTITIES C (1) , F (2) AND C (2) 

For time t ~ 1/e, we expand the evolution operator 
in At and collect the contributions that arise at sec- 
ond power in (A/e). By making use of the expression 
Eq. (|A12|) for the perturbation in the interaction picture 
we can compute the qubit components of the second or- 
der contribution to the continuous POVM. We define 

O s {t',t") = exp(ise(t' - 1") - is^{t' 2 -t" 2 )) 

x (6a(t')e lS9t '\5a(t")e lS9t "), (Bl) 

and (u\/3) is the overlap between coherent states, and 

eHt',n = fe«(io*+fe«(i"), (B2) 

Ci 2) (*'»*") = -SxW(t>)+5x^(t"). (B3) 

( 2^ 

The first term Q if , t ) represents the complex displace- 
ment of the oscillator position due to the perturbation 
acting one time at t" < t (forward in time), and one 
time at — t' > —t (backward in time). The second term 

Q (r,r ) represents the displacement of the oscillator 
due to the perturbation acting two times at t" < t' < t. 
Between the two perturbations the system evolves freely 
for the time t" — t 1 and accumulates a phase that de- 
pends on the difference of the effective qubit-dependent 
frequencies. In the short time approximation t « 1/e 
such a phase can be neglected. Integrating the position 
degree of freedom over the subsets ?y(s')i we obtain 
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F^{s\t) ss =s's^- [ dt' [ dt"Re\o s (t',t 
4 Jo Jo 



erf 



'Sx(t)+^ 2 \t',t")\ , J6x(t)+(i 2 \t',t") 



erf 



. (B4) 



where s = —s. This expression has meaning only in the short time approximation. By setting (» 1/e, the correction 
F^ 2 \t) at second order to the discrete POVM is evaluated to be 



2 i-t r i! 



4 Jo 



F^(t) = - \ dt' dt" cos (e{t' - t") - ip (t' 2 - t" 2 )) e 



12 - i" 2 \\ P -h\ A \ 2 (t'-t") 2 



{erf[|A|(t + t' + t")] + erf[|A|(t - t' + 1")}} 



(B5) 

In an analogue way we calculate the elements of the first and second order contribution to the double measurement 
operator C. The off-diagonal matrix element of the first order contribution C' 1 ' is 



CW(f; t) = \ (T(t) - 1) F«(i' -t) + ^crf ^E!LJL\ J* d t"e ut "T{t") 
and the full expression of the diagonal matrix element of the second order contribution is 
nV)M = - f j dtVW t '- t '')-^( t ' 2 - i '' 2 »e-^( t '- t '') 2 erf[|A|(t + t / + + 
F( 1 )(t)e^t , e-^ , r(t)rft / )*e-^ At 'l 2 +l A l 2t ' t erf ( ^^1 

I 



(B6) 



- Im 



(B7) 
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